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We have considered linear kinetic theory including the electron spin properties in a magnetized 
plasma. The starting point is a mean field Vlasov-like equation, derived from a fully quantum 
mechanical treatment, where effects from the electron spin precession and the magnetic dipole 
force is taken into account. The general conductivity tensor is derived, including both the free 
current contribution, as well as the magnetization current associated with the spin contribution. 
We conclude the paper with an extensive discussion of the quantum-mechanical boundary where we 
list parameter conditions that must be satisfied for various quantum effects to be influential. 
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I. INTRODUCTION 



Recently there has been much interest in the properties of quantum plasmas, see e.g. Refs. JlHi2- The research 
has been motivated by applications to e.g. quantum wells |12j . plasmonics [9l. Il3j. spintronics |14| , astrophysics 1 1 511 
and ultra-cold plasmas [16| ■ Quantum plasma effects has also been measured in solid density target experiments [17| . 
Within fluid theory, the theoretical models applied cover effects such as particle dispersion (the Bohm-de Broglie 
potential), and Fermi pressure. The magnetic dipole force and the magnetization current associated with the electron 
spin has also been captured within fluid models, with a macroscopic spin density as an extra dependent variable Q- 
More accurate models applying kinetic theory has also been studied [a, EH EH El • I n the absence of spin effects, a 
kinetic evolution equation can be derived for the Wigner function For long scale- lengths this equation reduces to 
^ \ the classical Vlasov equation. 

Including the spin degrees of freedom, it was recently shown [6| that a Wigner transform in regular phase space 
and a Q-transform in spin space produces a physically attractive evolution equation for a scalar distribution function. 
If the magnetic moment of the particles is put to zero, the ordinary Wigner equation is immediately recovered. Here 
we will study the opposite limit, however. For wavelengths much longer than the characteristic de Broglie length of 
Q_i, the particles, the quantum effects associated with particle dispersive effects (de-localized wavefunctions) disappear, 
and the evolution equation is much simplified. The remaining quantum effects are all due to the electron spin, 
, explicitly due to terms being proportional to the magnetic moment in the dynamic equations, or indirectly due to 

■ the background distribution function obeying Fermi-Dirac statistic. The kinetic model obtained from the long-scale 
limit of a fully quantum mechanical treatment is a slight generalization of a semi-classical spin model presented in 
Ref. ficjl. where the spin vector complements the regular phase-space variables as an independent variable (til Ho}. In 
Ref. 1 101 it was also shown that new wave modes appear in a magnetized plasma due to the combined dynamics of the 
magnetic dipole force and the spin precession. Here we will extend the linear analysis made in that work in several 

■ directions: 

1. By considering the case of general wave propagation (i.e. an arbitrary direction of propagation and a general 
— , wave polarization). 

^ . 2. By not making any restrictions on the wave frequency domain. 

• i-H ■ 

3. By using a somewhat more accurate model, derived from a fully quantum mechanical treatment. 

'_ ' 

C$ 4. By presenting a more thorough discussion of the classical - quantum mechanical boundary. 

The theory turns out to produce a Hermitian structure for the general conductivity tensor, including the free and 
magnetization currents, where all non-Hermitian contribution are associated with the poles. Contrary to common 
quantum effects associated with particle dispersion (de- localized wave functions), the effects due to the electron spin 
does not necessarily vanish in a plasma of moderate density and temperature. A low temperature and a high density 
tend to make spin effects more important (similarly to the quantum effects described in e.g. Refs. P, 0|), but a strong 
magnetic field or wave particle resonances can modify this picture considerably. As will be discussed in some detail in 
Section III, we point out that there is a large number of dimensionless parameters that together captures the relative 
significance of spin effects in different plasma regimes. Furthermore, in order to demonstrate the usefulness of the 
general conductivity tensor, in Section IV we present an analysis of the transition from wave damping to instability, 
when free energy is added to the the spin-degrees of freedom in the background distribution function. Specifically, it 
is shown that a background distribution function with equal population of the spin-up and spin-down states states 
leads to an instability. 
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II. THE CONDUCTIVITY TENSOR 



For long spatial scale-length the evolution of a (quasi) distribution function /(r, v,s, t) for electrons is described by 
the Vlasov-like equation [f| 
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where q e = — e is the charge of an electron, /i e = —(g/2)eh/2m e is the electron magnetic moment (note that the sign 
is included) and g = 2.002319 is the electron spin g factor, s is the unit spin vector, h is Planck's constant divided by 
27r, and m e is the electron mass. The distribution function / is normalized such that J fcrsdv = n with n being the 
number density. As for spin variables it is convenient to use ordinary spherical coordinates, i.e. cPs = siii9 s d6 s dip s . 

In the derivation of Eq. ([1} found in Ref. @ , trigonometric operators are expanded in powers of h. Consequently, 
Eq. ([IJ is the long scale limit of a more complete quantum description, where long scale here means that length 
scales are larger than the characteristic de Broglie wavelength of the electrons. Upon disregarding spin effects, the 
opposite limit of the full quantum theory, i.e. for short length scales, is the well-known evolution equation for the 
Wigner function [lj. The spin term proportional to s • B in Eq. JT]) is due to the magnetic dipole force, and the term 
proportional tosxB owes to the spin precession. Both these terms are also found in the semi-classical version of Eq. 
([T|), see Ref. The spin term proportional to B ■ Vg, however, is a fully quantum mechanical effect. Essentially 

this term can be viewed as a modification of the magnetic dipole force that occurs due to the spread out nature of 
the spin probability distribution. 

By solving Eq. ([T]) we may find an expression for the distribution function / from which we can construct the free 
current, Jf, and the magnetization current density, 3m, according to @ 

J = Jf + 3m 
= J/ + V x x M 



q e I vfd 2 sd 3 v + V x x I fj, e 



3sfd 2 sd 3 v 



(2) 



In this section we only consider the electron contribution to the current density. The ion contribution may typically 
be treated classically (because ion spin effects are suppressed compared to electron spin effects due the heavy ion 
mass), and its contribution may be found in standard literature, see e.g. Ref. (20| . The derivation given below may, 
however, easily be generalized to also apply for arbitrary spin-1/2 particles through the appropriate substitution of 
charge, mass and magnetic moment. 

Equations ([l) and ([2]) gives, together with the Maxwell's equations, a closed system of equations that can be solved 
to describe the dynamics of a plasma where quantum effects of the electron spin are captured. Below we will derive 
the conductivity tensor <Jij, defined as Ji — tiijEj, for such a system (2l| . With the conductivity tensor known, it is 
straightforward to construct a dispersion matrix and find the dispersion relations for arbitrary wave modes. Next, we 
linearize the kinetic Eq. ([T]) according to / = fa + fx and B = Bo + Bi, where the subscript denotes an unperturbed 
quantity and the subscript 1 denotes a perturbation, and we take Bo = BqZ. Before proceeding, let us point out a 
few quantum effects that may be contained already in the unperturbed distribution function: 

1. Fermi-Dirac statistics: This effect is well-known. Here we just note that for a plasma of moderate density 
and temperature this effect may be suppressed, as a large (negative) chemical potential (that applies for 
turns the thermodynamic equilibrium distribution into the classical Maxwellian. 



i^ 3 /m e k B T < 1) 



Landau- quantization: The quantization of perpendicular energy states becomes important in the regime of very 
strong magnetic fields, or very low temperatures, when h \uj ce \ /ksT 1, where ks is Boltzmann's constant and 
uj ce = —eBo/m e is the electron cyclotron frequency. 

Spin-splitting: The two spin states, up- and down relative to the magnetic field, clearly have different probability 
distributions in spin space. As a result, the general time-independent distribution function can be written as 
fa = /0+ + /0- with fo± — (l/47r).Fo±(v)(l±cos0 s ), where for a time-independent distribution function Fo± can 
be arbitrary functions of (v±, v z ), and Fg± is normalized such that J Fa±d v — n$± with hq± being the number- 
densities of the spin up/down states respectively. The positive spin state here means that the spin points in the 
direction parallel to the magnetic field, which means that the magnetic moment points in the opposite direction. 
Note that with this definition, the lower energy state is the spin state with negative index, i.e. no- > uq + in 
case the background distribution fo describes a thermodynamic equilibrium. 
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For a full quantum mechanical expression of the thermodynamic equilibrium distribution, see Ref. 0. However, in 
the present manuscript we will focus on the effects due to the dynamical equations, rather than the initial conditions 
briefly discussed above. Still, we note that since the spin-distribution has no classical limit point 3 above cannot be 
ignored. Since the present collision-free model has no mechanism for spin-flips to occur, we also note that the number 
of particles with spin-up N + = J fo+cPscPvcPr and the number of particles with spin down N- = J fo-d 2 sd 3 vd 3 r 
are conserved quantities. Furthermore, in thermodynamic equilibrium the spin-up and spin-down number densities 
no± = f fo±d 2 sd 3 v are related by («o+ — n Q-)/{ n o+ + ^o-) — tanh(/i e _Bo/fcs7 1 ). Since the evolution equation is the 
same, independent of the unperturbed spin state, it is convenient to wait with the split of /q into fo± to the final 
stages of the calculation. 

After linearization, Eq. (TTJ) is written as 



+ v ■ V X A + ^(vxB„). V v /i + %(s x B ) ■ V,/i 
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V v /o 



2^ £ 
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Next we make a plane wave ansatz of the perturbed parameters according to fi — j\ exp[i(k • x — uit)], etc. Without 
loss of generality, we define the wavevector as k = fc^x + k z z. We also choose to express the velocity in cylindrical 
coordinates (v±, ip v ,v z ) such that d 3 v = v±dv±dtp v dv z , and expand /i in eigenfunctions to the operator of the right 
hand side 



A = ^2 ^2 9ab(v±,v z ,d s )ip a (<p v ,v_ L )-^==exp(-ib<p s ), 
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2tt 
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where 



i>a(<Pv,v±) = 



/27T 



exp[-i(a(p v - kj_v± smtp v /uj ce )] 
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k±v± 



exp[i(l - a)ip v ] 



(5) 



where Ji{x) is a Bessel function of the first kind. We may then note the following simplifying relations; 

Qe , ^ tj x „ , dfi 
— (vxBoj-Vv/i = -w ce - — , 
m e ocp v 



2^ 



(s x B ) • Vg/i = -u c 



'dips' 



(6) 



where uj cg = 2^ p Bq/H is the spin precession frequency. Moreover, for simplicity we assume an isotropic distribution 
function on the form /o = /o(v 2 , 6 S )- For this case we may drop the term (v x Bi) • V v /o- We may also take advantage 
of the relations (df /dvj_)/2v±^ = (df /dv z )/2v z = df /dv 2 . 

Using the eigenfunction expansion of f\ (Eq. Q) in the linearized Vlasov equation ([3]), multiplying the resulting 
equation with ip*e iblps /-\/2w (where the star denotes complex conjugate) and integrating over <p v and tp s , we find the 
equation 



i(w - k z v z - au ce - buj cg )g ab = I a b{v±,v s ,(p s ), 



(7a) 



where 



lab = 



2-7T r 2ir 



dtp v dip s 



— E + ^V x (s • Bx + Bi ■ Vi) ) ■ V v /o + ^(s x B x ) ■ V,/, 

nip rrip \ ) h 



ipl—}==exp(ib(p s ). 



(7b) 

From Eq. ([Taj) we can then solve g a b and thereby construct /i in terms of Jo- It may be useful to note that the vector 
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products in (|7b[) arc; 

E • V v /o = 2 (E x v± cos<f v + E y v± smcp v + E z v z ) 
V x (s • Bi) • V v / ( 
V x (Bi • V s ) ■ V v /o 
(sxB^-Vg/o 



TO = 



dfo 



2i (k z v z + k±v± cos cp v ) (B x sin 9 S cos tp s + B v sin# s sin w s + B z cos S ) — 

Or 1 



2i (k z v z + k±v± cos <£>„) (i^ cos 9 S cos 99 s + cos 9 S sin ^ a — B z sin S ) 



(B x sin <p s - B y cos 



9/0 



(8b) 
(8c) 
(8d) 



With the relations ([5]) and (j8aj) - (|8d[) . the integration in j7b| can be performed in a tedious but straight forward manner 
(some useful integrals are found in appendix |A1 Eqs. (|Ala|) - (|Alc[) '), and through Eq. Q the resulting expression for 
/1 becomes 
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Here, the argument of the Bessel functions is k±v±/oj ce . In obtaining Eq. ([9]) we have used VxE = — <9B/9i to relate 
the magnetic field to the electric field, as well as the Bessel function identities J a +\{x) + J a -i{x) — 2aJ a (x)/x and 
Ja-x{x) - J a+1 {x) = 2Jl(x). 

With fi known in terms of fo, we proceed by splitting / into its spin states, f = fo+ + fo, and construct the free 
current 



Jf = ^2 q e J vfi v d 2 sd 3 v 
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as well as the magnetization current 
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where we have performed the tp v and (p s integration and immediately reintroduced the full integration element for 
convenience of a compact notation. We have also included the denominators in Q in A v by defining A v = A u /(uj — 
k z v z — auj ce + Lu cg ), and similarly for B v and C v . 

Before we construct the conductivity tensor, we note that we can simplify things yet further by writing fo± = 
(l/4-7r)(l ± cos9 s )Fq±(v 2 ). We may then perform the 9 S integrations in Eqs. (jTUJ) and (fTTj) (a useful integral is found 
in Appendix [Al Eq. (|A2|> ) ■ and construct the conductivity tensor which may be written as 
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where 



y-(ci) = o <ll dF 0y 
<»« m e dv 2 



is the classical contribution, and the spin contributions are 
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Given the conductivity tensor fj 12[) . the general dispersion relation is obtained in the same way as in the classical case, 
i.e. detDy = 0, with Dij = Sy — kikjC 2 /ui 2 ~ ia^/sou. Analogously to the classical case the conductivity tensor has a 
Hermitian structure (ov, = — <r|j) if the pole-contributions (associated with the denominators (ui — k z v z — acu ce — bui C g)) 
are dropped. It should be noted the Hermitian structure of the conductivity tensor does not follow trivially from Eqs. 
(TTU|) and (TTTj) . A detailed discussion on this is found in appendix [B] 

In appendix [Cl we derive the short Larmor radius limit of the conductivity tensor (|12p . i.e. the limit in which the 
argument of the Bessel functions are small, k±v±/uj ce <C 1. 



III. THE CLASSICAL - QUANTUM MECHANICAL BOUNDARY 

In addition to the spin phenomena studied here, using Eq. ([I}, kinetic quantum plasma phenomena associated 
with non- localized wave functions can be studied using the Wigner function, see e.g. Ref. Those quantum effects 
prove to be important for wavelengths shorter than the thermal de Broglie wavelength. However, collective effects 
are typically negligible for such short wavelengths, unless the parameter d± = huip/ksT^ approaches unity, where 
u}p = {riQe 2 / EQirif,) 1 / 2 is the plasma frequency. Here T m is the largest quantity of the Fermi temperature Tp and the 
thermodynamic temperature T. Thus it is safe to say that the quantum effects that are not associated with the spin 
requires either a low temperature and/or a high density. In addition to quantum effects associated with de-localized 
wave functions, there are those associated with the Fermi pressure. Naturally the relative importance of this effect 
is captured by the parameter 0I2 — Tp/T — (37r 2 ) 2 / 3 ft 2 n 2 / 3 /2m e fcsT, which also requires a low temperature and/or 
a high density plasma, although the scaling differs somewhat from that of d\. In linear theory, given a classical 
kinetic dispersion relation for a thermodynamic equilibrium system, the effects of the Fermi pressure is relatively 
trivial, since it can be deduced simply by replacing an unperturbed Maxwellian distribution with the Fermi-Dirac 
distribution. Estimating when the electron spin effects is of importance is somewhat more complicated, as the relevant 
dimensionless parameter now also involves the magnetic field strength and the wave-number and wave-frequency. As 
will be verified below, from the expression for the conductivity tensor components (|f 21) . we point out a set of important 
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parameters determining the relative strength of spin effects. Throughout this section we will assume the presence of 
ions described by a classical contribution to the conductivity tensor, see e.g. Ref. [20fl and electrons described by Eq. 

• The Zeeman energy over thermal energy, c?3 = fi e Bo/ksT = hui ce /kBT. This parameters enters for several 
reasons. The simplest one is that it determines the ratio between the spin-up and -down populations (i.e. Fo+ 
versus Fq- in (|12j) ) in thermodynamic equilibrium. 

• di = hk 2 /m e uj. This parameter is important since it determines the relative strength of the magnetic dipole 
force as compared to the force from the electric field. For systems that allow short-wavelength low-frequency 
modes this is of particular significance. 

• e?5 = hw ce /mic\, where ca = {Bq / Honiino) 1 ^ 2 is the Alfven velocity. This parameter can appear for several 
reasons. A relatively simple example is that it gives the relative strength of the magnetic dipole force to the 
Lorentz force, evaluated for ion-cyclotron waves. 

• e?g = hujp/m e c 2 . This parameter describes the relative magnetic permeability due to the spin, for a given spin 
state. When it approaches unity the density is so high that a relativistic treatment is called for. Thus the 
regime d@ > 1 is strictly not allowed by our treatment. However, as we will see below, the parameter do still 
have a certain significance, also for our non-relativistic case. 

Although the parameters d\ and d<i has not been directly associated with the electron spin, it should be noted 
that these two quantities can also be involved when determining the relative importance of spin for a specific wave- 
mode. In some cases it is a combination of the basic parameters d\ — do that gives the ratio of the spin and classical 
contribution. In order to illustrate the situation further, we will briefly discuss some specific examples. 

Magnetosonic waves: For usual ideal MHD waves there are the shear Alfven wave with E x ^ 0, and the fast 
and slow magnetosonic waves with E y>z ^ 0. The latter wave mode is determined by the components a yz (= —v* zy ) , 
a yy and a zz . Assuming that the plasma is strongly magnetized or has a low temperature, we consider the case when 
c?3 approaches unity, but all other quantum parameters are small. Considering the MHD-regime it is then found that 
the ratio of the spin contribution to a yz (from (Y^j)) to the classical contribution (from (ly )) is o^f 1 /a y °P ~ d§. 
For a low-beta plasma this mainly affects the slow magnetosonic mode, but for a high beta plasma both the fast and 
slow magnetosonic mode is much affected in this regime. This example mainly applies to strongly magnetized plasmas 
that occur in astrophysics, e.g. pulsars and magnetars. If experiments on ultra-cold plasmas [16| were extended to 
include strongly magnetized plasmas, this could also be a regime of relevance. 

Ion-cyclotron waves: When ordinary MHD waves (e.g. shear Alfven waves and compressional Alfven waves) have 
their wavelengths shortened to k z ~ io c i / ca and k^ ~ u c i / ca where cu c i is the ion-cyclotron frequency, the frequencies 
are increased up to lj ~ uid, and the waves become dispersive. For arbitrary directions of propagation, all components 
of the conductivity tensor are of relevance, although the classical expression for <r xz (and a zx ) is usually possible 
to put to zero. For definiteness we will make a comparison of a xx , but other components would work as well. We 
consider a regime of relatively high density, but with a magnetic field that is not particularly strong. For example, 
no = 10 21 cm -3 and Bo = 10~ 4 T. We note that this regime means that c\ <C c 2 and for a moderate temperature 
(such that the Bessel functions can be Taylor expanded to lowest order) the only large quantum parameter is d$ (and 
possibly di). Evaluating the classical part of a xx (that comes from the ions) and the spin part of the same component 
(that comes from the terms proportional to X^ p ^(or with a = ±1 such that the relevant denominators in (fT2)) 

becomes uj — k z v z ± (uj ce — uJ C g), we find the ratio <Txx^ /&xx = d 2 , where we have used that |w ce — uj cg \ is of the order 
of the ion-cyclotron frequency uj C i. In general all of the different ion-cyclotron modes are much affected when d$ — > 1 

Whistler waves: Next we focus on the case of whistler waves, lu < |w ce |, where the electron conductivity is 
dominant also for the classical contribution. Furthermore we let fcj_fth/ \w ce \ 1: such that the Bessel functions 
can be expanded, and we note that we are considering the regime of short perpendicular wavelengths, k± < \uj pe \ /c. 
The most relevant components here are a xx , a yy as well as a xy . The dominant spin contribution comes from the a xy 

component of ^ sp '' , which should be compared to the a xy component of Y i ^ . Using the thermodynamic equilibrium 

relation (rio+ — n.o-)/(^o+ + n>0— ) = tanh(/i e i?o/fcBT) rs fi e Bo/kBT, we then find the ratio a x s y ^ /cr x c y — d^d^, and 
similar comparisons can be made for other conductivity components in the whistler regime. Note that for a non- 
relativistic Alfven velocity we can simplify d^d^ — de(hio p /kBT). Since it is the thermodynamic temperature that 
enters here (which determines the relative population of the two spin states) , rather than the Fermi temperature, we 
can very well be in the non-relativistic regime de 1 simultaneously as dQ(hio p /kBT) ~ 1. For example, a plasma 
with density no — 10 30 cm 3 and T = 10 5 K is nonrelativistic to a good approximation (i.e. do <C 1) at the same time 
as d^ds ~ 1. 
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The above brief examples illustrates some possibilities to get conductivity components due to the spin that are 
comparable to the classical ones in magnitude. Note that it is hard to give simple general guidelines when spin effects 
are important. For example, for a given temperature we may need a sufficiently strong magnetic field for spin to be 
important. We saw this in the example with magnetosonic waves, where the parameter ^3 should be of order unity. 
On the other hand, for a given density, and for certain wave modes, we may need a sufficiently weak magnetic field 
for spin to be important, as we saw in the example with ion-cyclotron waves involving the parameter d^. 

Furthermore, sometimes the spin terms are small compared to the classical terms, but still give a significant 
contribution. The reason is that a number of small but qualitatively new effects from the spin can be added to the 
dispersion relation. Thus, in order to get a somewhat deeper understanding of the classical - quantum mechanical 
transition, than mere estimates of magnitude can give, we will present a specific calculation in more detail. In the 
example below, we follow the assumptions made in Ref. [ll|. We then focus on the regime fcx^tft/ l w ce| <S 1) such that 
the Bessel functions can be expanded and consider c?3 <C 1. As pointed out previously, the latter strong inequality 
condition holds for most plasmas, except for some strongly magnetized astrophysical objects [15j . Furthermore, to 
be specific, we consider waves with a polarization E = E y y and B = B x 5i + B z z = (k z E v /uj)Sc — (k±E y /u>)z. For 
this polarization to be possible we must have cr xy , o~ xz <j yyi which can be verified a posteriori. We keep the 
standard classical terms (see e.g. Ref. |20j) up to order uj~ 2 in an l/w ce -expansion (for both ions and electrons), where 
the dimensionless expansion parameters are considered to be io/uj ce , k z Vth/u ce and k_\_v t h/u) ce . The spin effects are 
assumed to be smaller, due to H \ui ce \ jkuT <C 1, and we also consider the case uj / \lu cc — Lu cg \ < 1. Accordingly we 
only include spin terms to zero:th order in the l/w ce -expansion. With these prerequisites the only spin effect that 
survives is the z-component of the magnetic dipole force, and its corresponding modification of the y— component of 
the magnetization current. For the given polarization the dispersion relation then reads 



lu 2 - k 2 c 2 - —o-yy = (13) 

£0 



where a yy is given by 



zeoW _ +zeo ___ J _____ dVz (14) 



CI 



Asxrnil J (w~k z v z ) 

where all integrations except that over dv z has been carried out, and the normalization of the rescaled distribution 
function Fq is given by J F$dv z = 1. In Eq. (|14| the first term of the right hand side is the classical contribution from 
the ion free current, whereas the second term is due to the electron magnetization current. The dispersion relation 
(fT3")) now immediately reduces to 



k 2 c 2 



■uj 2 ./oj 2 

VII ■ 



(15) 



. 2 2 H 2 k 2 UJ 2 e f Fq 

sin a cos a— — rHr- / ^dv z 

J {LO~k z V z f 

where sin 2 a = k\jk 2 . Firstly, we note that for ordinary densities in space and laboratory plasmas, the influence 
of the real value of uj is small and thus we can omit the spin term as well as the first term of Eq. (|13p and write 
the standard dispersion relation for compressional Alfven waves lu 2 = k 2 c\. where we have assumed ca *C c (or 
Upi/t^ci ^1)- Furthermore, for Re (a;) to be significantly changed by the spin effects, we need high densities and/or 
low temperatures such that h 2 uj 2 e /m 2 c 2 v 2 approaches unity. However, increasing the density to fulfill this, the 
minimum velocity spread will also increase, and when this parameter approaches unity we have also reached the 
regime of a relativistic Fermi velocity. Nevertheless, one can see that the spin has a certain significance even if the 
dimensionless parameter of this example is much smaller than unity. Firstly, the spin term makes the dispersion 
relation slightly anisotropic, since the group velocity is not parallel to the wave vector anymore. Secondly, in the 
given approximation the group dispersion comes solely from the spin effect. Finally, provided the omitted cyclotron 
resonances lies far out in the thermal tail, the imaginary contribution is determined by the spin term, and is given by 

■ ' ' x 1 2 ^sin 2 ah 2 ^ /(Rew) 2 
1 + ijj 2 Juj 2 t J 2m 2 c 2 " 2 6XP ' h2 " 2 



7 = Im(w) = - ; , 7-2 .... ■. , .■ ex P nor ( 16 ) 



when the influence of the spin term is small. Thus we can conclude that there are several possibilities for spin to 
be of significance, even if all dimensionless parameters determining the magnitude of the spin contribution are much 
smaller than unity. 

As a final comment, let us make a few comments regarding the terms in Eq. (|12p with denominators uj — k z v s ± 
{uj ce ± uj cg ). As pointed out in Ref. [loj], these terms may give raise to wave modes with frequencies 
which is a completely new effect due to spin, that can survive even when all dimensionless parameters d\ — d§ are 
small. In the next section we will demonstrate from the general result, Eq. (1121) . that the presence of the denominators 
uj — k z v s ± {uj ce ± ui cg ) has further consequences for the transition from wave damping to instability. 



IV. SPIN INSTABILITY - AN EXAMPLE 



In order to illustrate the usefulness of (|12|) , we will evaluate the imaginary contribution associated with the denom- 
inators u) — k z v z ± (uj C g — ui c ). To compute such terms explicitly, we must pick a specific form of the unperturbed 
distribution function. As we here have not been interested in classical instabilities, that can be induced by a velocity- 
space non-equilibrium, we have assumed the distribution function to be isotropic, i.e. fo = fo(v 2 ). However, as we will 
see below, this does not remove all possibilities for free energy to be present in the background distribution function. 
Specifically, in case the number of particles in the two spin states does not correspond to thermodynamic equilibrium, 
we may find that certain wave modes may become unstable. To illustrate this idea we let the unperturbed distribution 
function Fq± be of the form 



3/2 

F ± = n [ "~ c ) exp 



m e v 2 



2fcs7kin 



G±(T sp ) 



where G±(T sp ) is proportional to the number of particles in spin- up and spin-down states respectively and can be 
written as 

q r T )= cxp(vfi e B /k B T sp ) 

sp eyLp(fi e B /k B T sp ) + exp(-fi e B a /k B T sp ) 

where we have introduced the parameter T sp that can be interpreted as a spin temperature, in addition to the kinetic 
temperature T^m of the ordinary Maxwellian velocity dependence. If T sp = Tki n = T we get a thermodynamic 
equilibrium distribution with the common temperature T, and if we let T sp — ¥ oo, we get the same number of particles 
in the two spin states. We note that for most plasmas, the energy difference between the high-energy and low-energy 
spin state is small ( i.e. typically fi e Bo/k B T <C 1), but nevertheless even a small energy energy difference can be 
crucial, as we will see below. 

Next we assume that the plasma parameters (Bo, no, Tk; n ) correspond to a classical regime, such that the real part 
of the frequency co r (k±, k z ) is given by the classical terms in (|12[) to a good approximation. Furthermore, we assume 
that all resonant electron velocities v z = ±(w — nu c )/k z corresponding to the classical terms lies very far out in the 
thermal tail. The dominant imaginary contribution to the dispersion relation may then come from the spin terms 
x\j P ^ and Z^j P ^ with resonant particle velocities v z = (uj ± \Auj c \)/k z (corresponding to a = ±1 in the conductivity 
tensor (H2)) with |Aw c | = \oj cg — u) ce \. Furthermore, letting the Larmor radius be smaller than fc^ 1 (such that the 
Bessel functions can be Taylor expanded) , we can compute the pole contributions from X>j P ' and Z-j with the help 
of well known properties of the plasma dispersion function; ImZ(£) = in 1 / 2 exp(— £ 2 ). Since the lower value of the 
resonant velocity occurs for v z — (uj — \Auj c \)/k z we concentrate on the contribution from x!-j (the contribution 
from z\° p) at v z = (u) + \Aw c \)/k z can be computed analogously). The result is then 



/27rfc B T kin yL 2 e k\ ( ( fi e B \ , fko ce , 

Im(<7jj) = n Q \ — 5— -tanh - — — — + — — — — exp 



Ahk z w 2 uj \ \k B T sp J 2k B T k 



-m e (lo - |Au; c |) 2 
2k 2 k B Tun 



= Aij (T kin , n , B ) (tanh (|aT) - a) (17) 

where we have noted that "Yli V G v = 1 and vG v = tanh(fj, e Bo/k B T sp ), and we have defined a = 2fcfS- an< ^ 
T = 7kin/7sp. In Fig. 1 we plot Im(<7ij) as a function of T for some different values of a. For the case of thermodynamic 
equilibrium, T sp = Tkim the fact that g/2 > 1 ensures that Im(aij) > (at least for not too strong magnetic fields, 
i.e. small values of a). In this case the contribution to the conductivity tensor from the poles gives raise to a small 
damping rate. This applies independently of the polarization of the mode under consideration. However, we need 
only perturb the spin temperature a little from equilibrium (T sp > T^ n ) for the system to become unstable. The weak 
damping then turns to a weak growth rate of all modes where the resonant velocities v z — (to — Aw c )/fc z does not 
fall too far out in the tail of the distribution. We stress that this may apply to all sorts of ion-cyclotron modes, as 
A.u> c is of the same order as the ion-cyclotron frequency. It is not hard to imagine situations where the plasma state 
deviates from thermodynamic equilibrium, in which case the sign of the pole contribution may change. For example, 
if particles are flowing from a region of weak magnetic field to a region of strong magnetic field, the inflowing particles 
will be more evenly distributed between spin-up and spin-down states than that of the local thermodynamic state. 
Hence the spin-distribution will be more evenly distributed than for T sp — T^ ln , i.e. we will get T sp > T^ n . The 
same situation will also occur if the external magnetic field is increased slowly, but at a rate slightly faster than the 
spin-states can relax to thermodynamic equilibrium. 
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FIG. 1: (Color online) The imaginary contribution of X\° p to cry plotted as a function of T sp /Tki n for different values of a. A 
positive value corresponds to damping while a negative value gives rise to a instability. 



It may be noted that if a is sufficiently large (corresponding to a strongly magnetized systems), the result (fl7|) 
predicts that the system will be unstable already at thermal equilibrium. However, as a increases, the effect of 
Landau quantization will become important (which affects the perpendicular kinetic energy distribution, see e.g. |6j|), 
and hence the background distribution function will no longer be isotropic. The validity of Eq. (JTTJ) is therefore 
restricted to small values of a. Finally, let us make an estimate for the growth rate for the case of ion-cyclotron Alfven 
waves with uj ~ uj C i ~ k z CA ~ k±CA- We let the the temperature ratio be well above the threshold for instability (i.e. 
T S p/Tki n ~ 2), and let a -C 1 (such that the isotropic distribution function is justified). The normalized growth rate 
7/cj found from [T7] is then of the order 

7 fi 2 wL , , 

*V 18 

u> m e miC A 

if the let the resonance approach the bulk of the distribution in order to find the maximum groth rate. 



V. CONCLUSION 



In the present paper we have studied the regime of long spatial scales in an otherwise fully quantum mechanical 
kinetic model. The assumption of long spatial scales (compared to the characteristic de Broglie wavelength) makes 
the model reminiscent of semi-classical theory. In particular, the distribution function behaves as if particles are 
effectively localized spatially, and unlike the Wigner function the distribution function is always positive in our case. 
Still, the impossibility for the spin vector of individual particles to have a localized probability distribution in spin 
space keeps certain quantum features in the evolution equation JT]). The linearized theory is solved in a magnetized 
plasma for a homogeneous background, which leads to the general conductivity tensor, Eq. (|12p . that includes both 
the contribution from the free current density as well as that from the magnetization current due to the spin. This is 
the main result of the paper. The main restriction in applying Eq. (|12| is that the velocity distribution is assumed 
to be isotropic. Specifically we have used f = ./o(« 2 ,#s) = (1/4tt) [F 0+ (v 2 )(l +cos9 s )+ T _(w 2 )(1 — cos# s )]. Note 
that, by contrast, the spin distribution of the unperturbed distribution function is the most general time-independent 
solution for a constant magnetic field background. 

In section III we have discussed the quantum-classical boundary, which is considerably more complicated when 
spin effects are included. In particular, as deduced from the dimensionless parameter c?3 — de, spin effects may in 
certain cases remain also in a plasma of modest temperature and density. Furthermore, effects due to the new types 
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of wave-particle resonances may be of significance independent of the dimensioneless parameters. In section IV, we 
have found that such resonances is one of the key ingredients in a new type of instability, somewhat reminiscent of 
the Weibel-instability. As is wellknown (see e.g. [20|), in the case of Weibel-instabilites the instabilities arise from 
a difference in perpendicular and parallell kinetic temperatures. Here, the deviation in the spin temperature from 
a common (isotropic) kinetic temperature is the source of the instability. Furthermore we note that a very small 
deviation of Tki n /T sp from unity is sufficient to drive an instability. 

Although the present result (jT2J) is relatively general, there is still several interesting extensions that could be made; 

1. Generalization to an arbitrary (non-isotropic) background velocity distribution. 

2. Inclusion of the spin-orbit coupling. 

3. Using the full evolution of Ref. [fj, including the short-scale physics (i.e. shorter than the characteristic de 
Broglie wavelength). 

4. Inclusion of collisional effects. 

Furthermore, a thorough evaluation of Q12p , that must be done numerically, could possibly reveal new and interesting 
possibilities, involving e.g. new wave modes. 

Appendix A: Some useful relations 

The integrals (|AlaP - (|Alcp are very useful when performing the ip s and ip v integrations in Eqs. (|7b[) . (fTT)|) and (fTTj) . 



2jv 







2ji 



, , 7r for n = ± 1 / a 1 % 

dtp exp{intp) cos(^) = <j Q foj . n ^ ±l) ( Ala ) 



, i ±wr for n = ±l / 4 n\ 

dtpexp{imp)sm(tp) = <j Q for n ^ ±l) ( Alb ) 



2 7T 



, , 2tt for n = , . .. , 

c^exp^) = <j Q for n _£ (Ale) 



When performing the 6* s integration in going from Eqs. (|10I) and (fTTIl to (fT2|) . it is useful to note that integrals on the 
form 



d6 (cos 6>)" (sin^)" 1 (A2) 
vanish for all odd integer values of n independent of the integer value of m. 



Appendix B: The Hermitian structure of the conductivity tensor 

The Hermitian structure of the conductivity tensor (er^ = — <j\ ) is not transparent from Eqs. (fTU|) and (fTTj) . The 
conductivity tensor will for instance not be recognized to have a Hermitian structure until after the (^-integration is 
performed. Below we list some further identities that have been used in obtaining cr,j on the form written in Eq. (1121) . 

To match the terms a xy with <j yx , it is found that an integral of the form 



aJl 9fo v ± dv ± (Bl) 



must match the integral 



^— In U! — k z V, — ClljJce 8v 2 



a {auj ce + k z v z )Jl df 
/ I i -7r^ v ±dv_L 

z — ' In LO UJ — K Z V Z — aUlr.p. OV^ 
a— — oo 



a— — oo 



/ , , «J„ + o^wxdwx- (B2) 
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For the terms to match, the following integral must vanish; 



i=— oo u a— — oo u 

^ /»CO OO 

= ~o / ^ aJa(Ja-l - Ja+l)fodV± 



a— — oo 

/•OO oo 



^ /'CO ^ 

•'° n=l 



2 

= (B3) 

where we have integrated by parts, summed over all integers a and used the relation J- n = (— l) n ^ n . Thus, since 
the integral (|B3|) indeed vanish, we have a xy = —c* as expected. 

In a similar manner, we find that the following integral must vanish for a yz to match a zy ; 

v z ^dv z = Q. (B4) 

This is trivially satisfied under the assumption fo(v 2 ) since this is a product of an odd and an even function integrated 
from — oo to oo. 



Appendix C: The short Larmor radius limit 



Below we consider the short Larmor radius limit of the conductivity tensor (|12p . i.e. the case where the Bessel 
argument is small k±v±/ \co ce \ <C 1. Since we have assumed an isotropic distribution for the plasma (fo(v 2 )), we note 
that 

dF _ 1 dF 
dv 2 2v z dv z 

With this noted, we can integrate by parts with respect to v z to eliminate the derivatives of Fq in the conductivity 
tensor (fT2]) . We encounter integrals on the form 



vj dFp dy = 1 
ft — k z v z dv 2 z 2 



dF 



Q-k z v z dv z dVz 
'{n-l)v^- 2 



k z vl 



Q - k z v z (Q - k z v z y 



dv z 



(CI) 



where n may take the integer values n = 0, 1, 2. Here Q = oj — aoj ce — bu cg with the integer values of a and b not being 
specified. Next, we expand the Bessel functions in v±kj_/uj ce and perform the summation over all integer values of a. 
Terms proportional to kj_vj_/uj ce are considered small and are dropped. The conductivity tensor then becomes: 



E [ 



.(*P) 



*{v)ij (v)ij 



(C2) 



v=+, 



where 



(cl) 
V(v)ij 



m e 



(u — k z v z ) 2 



/ \ 


V i J 

is the classical contribution, and the spin contributions are 



1 \ u}-2kzV z -au C e v± 
4 / j (uj — k z v z —au) ce )' 2 v 2 
a=±l 



y(«P) 



?0v 



(u — k z v z ) 2 



I -i a \ 

—a — i 
V / 

\ 



d 6 v 







Vo 



k\k\ ■ i v 2 ± k 2 ± (u:-2k z v z ) 



'- V I 



2 v 2 



vk i k 7 



-vk±_ k z 



/ 
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a=±l 



(ui — k z v z —auj ce ) 2 v z ^ 



/ 



V o 



-va 



ia—— + vi- — 



\ 



0/ 



d A v 



together with 



„(sp) i z (sp) 



£ 

b=±i 



F 



0// 



-vb 



1 ?J 



w - k z v z - buj cg 2m(w - fc^ - frw cff ) 2 



1 ft k\ v\ 



aui, 



bu, 



eg 



M^d 3 v 



a'fil 8 TOe Wce ^ ( W ~ kzVz ~ aUJce ~ buJc 

where 1 ' = My with My being defined as before, and M^- = M*^. 

We may note that we only encounter four kinds of unique integrals in the conductivity tensor (|C2[k 



F 



1 



Q - k z v z 
v 2 , 1 



d 6 v, 



;d 3 V, 



vl (ft - k z v z ) 2 
Fq^—— — t —d 3 v, 



F Q 



v z (Q - k z v z ) 2 

' ;d 3 V. 



(o - k z v z y 



For the special case when the cyclotron resonances are far out in the thermal tail, the conductivity tensor (|C2[) may 
be significantly simplified. This is, however, trivial to obtain from ()C2[) and will therefore not be pursued further here. 
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